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WASP-3b: a strongly-irradiated transiting gas-giant planet 



o 
o 

(N 
> 

o 



D. Pollacco 1 *, I. Skillen 7 , A. Collier Cameron 2 , B. Loeillet 11 , H.C. Stempels 2 , 
F. Bouchy 12 ' 13 , N.P. Gibson 1 , L. Hebb 2 , G. Hebrard 12 , Y.C. Joshi 1 , I. McDonald 5 , 
B. Smalley 5 , A.M.S. Smith 2 , R.A. Street 1 ' 15 , S. Udry 10 , R.G. West 3 , D.M. Wilson 5 , 
P.J. Wheatley 9 S. Aigrain 6 , C.R. Benn 7 , V.A. Bruce 2 , D.J. Christian 1 , W.I. Clarkson 4 ' 14 , 
B. Enoch 4 , A. Evans 5 , A. Fitzsimmons 1 , C.A. Haswell 4 , C. Hellier 5 , S. Hickey 7 ' 16 , 
S.T. Hodgkin 6 , K. Home 2 , M. Hrudkova 7 ' 17 , J. Irwin 6 , S.R. Kane 7 , F.P. Keenan 1 , 
T.A. Lister 2 ' 5 ' 15 , P. Maxted 5 , M. Mayor 10 , C. Moutou 11 , A.J. Norton 4 , J. P. Osborne 3 , 
N. Parley 4 , F. Pont 10 , D. Queloz 10 , R. Ryans 1 , and E. Simpson 1 



1 Astrophysics Research Centre, School of Mathematics & Physics, Queen's University, University Road, Belfast, BT7 INN, UK 
2 School of Physics and Astronomy, University of St Andrews, North Haugh, St Andrews, Fife KY16 9SS, UK 
3 Department of Physics and Astronomy, University of Leicester, Leicester, LEI 7RH, UK 
4 Department of Physics and Astronomy, The Open University, Milton Keynes, MK7 6AA, UK 
5 Astrophysics Group, Keele University, Staffordshire, ST5 5BG 
i 6 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CBS OB A, UK 
Oh 7 Isaac Newton Group of Telescopes, Apartado de Correos 321, E-38700 Santa Cruz de la Palma, Tenerife, Spain 

8 Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 3261 1-2055, USA 

9 Department of Physics, University of Warwick, Coventry CV4 7AL, UK 

10 Observatoire de Geneve, Universite de Geneve, 51 Ch. des Maillettes, 1290 Sauverny, Switzerland 
^ ■ 11 Laboratoire d'Astrophysique de Marseille, BP 8, 13376 Marseille Cedex 12, France 

12 Institut d'Astrophysique de Paris, CNRS (UMR 7095) - Universite Pierre & Marie Curie, 98 bla bvd. Arago, 75014 Paris, France 
13 Observatoire de Haute- Provence, 04870 St Michel V Observatoire, France 
T—( , 14 STScI, 3700 San Martin Drive, Baltimore, MD 21218, USA 

J> ■ 15 Las Cumbres Observatory, 6740 Cortona Dr. Suite 102, Santa Barbara, CA 93117, USA 
\J^) ' 16 Centre for Astrophysics, Science & Technology Research Institute, University of Hertfordshire, Hatfield, ALIO 9AB, UK 
i 17 Astronomical Institute, Charles University Prague, V Holesovickach 2, CZ-180 00 Praha, Czech Republic 



o 



c3 



O 



Accepted 1988 December 15. Received 1988 December 14; in original form 1988 October 11 



ABSTRACT 



k>( \ We report the discovery of WASP-3b, the third transiting exoplanet to be discov- 

ered by the WASP and SOPHIE collaboration. WASP-3b transits its host star USNO- 
B1.0 1256-0285133 every 1.846834 ± 0.000002 days. Our high precision radial-velocity 
measurements present a variation with amplitude characteristic of a planetary-mass 
companion and in-phase with the light-curve. Adaptive optics imaging shows no ev- 
idence for nearby stellar companions, and line-bisector analysis excludes faint, un- 
resolved binarity and stellar activity as the cause of the radial-velocity variations. 
We make a preliminary spectroscopic analysis of the host star finding it to have 
T e g = 6400±100K and log g — 4.25±0.05 which suggests it is most likely an unevolved 
main sequence star of spectral type F7-8V. Our simultaneous modelling of the transit 
photometry and reflex motion of the host leads us to derive a mass of 1.76^0 

and radius 1.31 j^' I^Rj for WASP-3b. The proximity and relative temperature of the 
host star suggests that WASP-3b is one of the hottest exoplanets known, and thus has 
the potential to place stringent constraints on exoplanet atmospheric models. 

Key words: methods: data analysis - stars: planetary systems - techniques: radial 
velocities - techniques: photometric 
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1 INTRODUCTION 



(200 



Since the discov ery 
ICharbonneau et al.l 
found to transit its host 
ther 22 transiting systems 



55 



iHenrv et al.l <|2000l ) and 

of the first exoplanet 
star, HD209458b, a fur- 
have been announced (see 
http://obswww.unige.ch/~pont/TRANSITS.htm). Tran- 
siting exoplanets are highly prized because the transit 
geometry constrains the orbital inclination, and this in turn 
allows their masses and radii to be determined directly. The 
mass-radius relation for exoplanets allows us to probe their 
internal structure, since it is these parameters which are 
compared with models of p l anetar y structure and evolution 
|Sato et al.1 120051 ; iGuillotl 120061 ). The limited numbers 
of transiting exoplanets studied so far show remarkable 
diversity in their physical parameters. For example, planets 
with masses M ~ lMj range in size from 0.8-1.5i?j for 
reasons that still elude us, alth ough several plausible ex- 
planations have bee n proposed (fee Arras fc Bildstenl 120061 : 
iBurrows et all 120071 ; iFortnev et al.1 120071 ; iGu et all 1200-3 ). 
The discovery of transiting planets in greater numbers 
will allow us to further explore the mass-radius plane, 
and thereby constrain theories of planetary formation, 
migration and evolution. 

The discovery of the firs t exoplanet, 51 Peg b, w as by 
the radial- velocity method (Ma yor fc Queloa |l995), and 
this technique is responsible for the discovery of the vast 
majority of the known exoplanetary systems, including 
HD209458b. Howev er, following the d i scovery of the tran- 
sits o f HD209458b (|Henrv et al.l l200Ct ICharbonneau et all 
l200d ) it was widely believed that the multiplex advantage of 
wide-field photometric imaging could lead to this technique 
becoming the dominant method for detecting exoplanets. 
While it is true to say that, at least initially, photomet- 
ric survey s have been slow to realise their expected detec- 
tion rates (|Hornell2003l ). recently this situation has begun to 
change with 14 new systems published in 2006-07 alone. This 
improved detection rate is largely due to the development 
of a better understanding of noise characteristics, especially 
the correlated no ise inherent in such photometric surveys 
jPont et alj|200d ). Transiting planets now comprise approx- 
imately 10% of the known exoplanets. Successful exoplanet 
photometric surv eys include Wide An gle Search for Planets 
(WASP) Project |Pollacco et alj2006l)7the Hungari an Auto- 
matic Telescope (HA T) Network jBakos et al.|[20o3 ). OGLE 
dUdalski et alJl2002h. the Tr a nsatlantic Exoplanet (T rES ) 
Survey l|Dunham et al.1 12004 lO'Donqvan et all 120061 ) and 
the XO group jMcCullough et alj|2006t ). TheWASP project 
has published t wo new systems in the last y ear WASP-lb 
and WASP-2b jCoUier Cameron et al.ll2007al ) and WASP- 
lb, in particular, has proved to be especially interesting, 
being well oversized compared to other planets for its mass. 

In this paper the WASP and SOPHIE collaboration 
announce the discovery of a new, relatively high mass, 
strongly-irradiated gas-giant exoplanet, WASP-3b. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 SuperWASP-N Photometry 

The photometric observations used in this study are from 
the inaugural 2004 SuperWASP-N observing season, which 
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Figure 1. Light Curves for 1SWASP J183431. 62+353941.4 
(WASP-3) obtained with SuperWASP-N (top panel), the IAC 80- 
cm telescope (V and I) and the Keele 60-cm telescope (R) as 
marked. All the data (apart from that from SuperWASP-N) 
was averaged in 300 second bins. The data was phased using 
the ephemeris derived in Section 13.2.11 T = 2454143.8504 and 
P = 1.846834 days. 



ran from April to September of that year; this data set also 
led to the discovery of WASP-lb and WASP-2b. Briefly, 
in 2004 the SuperWASP-N instrument comprised 4, and at 
times, 5 optical cameras, each consisting of a Canon 200 mm 
f/1.8 telephoto lens imaging onto a thermoelectrically- 
cooled, science-grade 2048 x 2048 CCD camera (manufac- 
tured by e2v Technologies PLC). In this system, the CCD's 
13.5 iim pixels project to an angular size of 14.2 arcseconds. 
For the entirety of the 2004 season, while robotic operation 
was being commissioned, the instrument performance was 
supervised with an observer always in attendance, although 
data acquisition was fully automated. Data were shipped to 
the UK on a weekly basis and reduced at the home institutes 
of the WASP Consortium using a dedicated, purpose-built 
pipeline, and the results ingested into the project database 
at the University of Leicester. The entire project infr astruc- 
ture is described in detail by iPollacco et all (|2006l ) along 
with the deployment of a further facility, WASP-S, at the 
South African Astronomical Observatory. SuperWASP-N 
now runs completely robotically, and data are transferred 
to the UK in near real-time over the Internet. 

Transit sea r ches were carried out on this dataset 
^Christian et al.l 120061: IClarkson et al.1 120071 : iLister et all 



l2007t IStreet et al.l I2OO7F" using the tec hniques out- 
lined by ICollier Cameron et all l|2006h . 1SWASP 
J183431. 62+353941.4, whic h we hence f orth denote as 
WASP-3, was highlighted bv lStreet et al.1 (|2007l ) as a high- 
priority candidate worthy of further study. Figure [1] (top 
panel) shows the original SuperWASP-N lightcurve, which 
comprises 3969 data points obtained over a 118 day period. 
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In the original SuperWASP-N photometry 17 transits were 
observed with >50% of a transit observed on 10 ocassions. 
These data led to an ephemeris of T' =2453139.1748 and 
P—l. 846800 which was used to arrange followup observa- 
tions. The transit here has a depth of 0.013 mag. and is 137 
minutes in duration. 



2.2 Higher precision photometric observations 

WASP-3 was observed with the IAC 80 cm telescope as part 
of the Canarian Observatories' International Time Pro- 
gramme for 2007. The imaging camera on this telescope has 
an e2v Technology PLC CCD of 2148 x 2148 pixels giving 
a scale of 0.33 arcseconds/pixel and a total field of view of 
10.6 arcminutes. Observations were taken during the transit 
of 2007 August 4, and consist of 327 images of 30 and 20 
seconds integration in the V and / bands respectively. This 
night was photometric but suffered from significant Saha- 
ran dust extinction, estimated to be ~ 0.4 mag on La Palma 
from the SuperWASP-N real-time pipeline. 

The images were bias subtracted with a stacked bias 
frame and flat-fielded with a stacked twilight flat field im- 
age obtained in both filters using individual flats gathered 
over the course of the run. After the instrumental signatures 
were removed, source detection and aperture photometry 
were performed on all sci ence frames using th e CASU cata- 
logue extraction software ijlrwin fc Lewidl200lh . We chose an 
aperture size matched to the typical seeing (5 pixels, 1.5") 
and selected 5 non-variable comparison stars in the field of 
WASP-3 to use in deriving the differential photometry. For 
each exposure, we summed the fluxes of the 5 comparison 
stars and divided by the flux of the target star to derive the 
differential magnitude of the target. The resulting V and I 
band lightcurves (Figure [TJ of WASP-3 have a precision of 
~ 4 millimag. 

Further observations of WASP-3 were made with the 
Keele University Observatory 60cm Thornton Reflector on 
2007 September 10. This telescope is equipped with a 765 
x 510 pixel Santa Barbara Instrument Group (SBIG) ST7 
CCD at the f/4.5 Newtonian focus, giving a 0.68 arcsec- 
ond/pixel resolution and a 8.63x 5.75 arcminute field of 
view. During most of the period the weather was photo- 
metric except post egress where some cloud appeared. Al- 
together 644x20 sec observations in the R band were ob- 
tained. After applying corrections for bias, dark current and 
flat fielding in the usual way, aperture photometry on two 
compariso ns were performed usin g the commercial software 
AIP4Win ijBerrv fe Burnellll2005r i. 

Tracking errors and spurious electronic noise mean that 
systematic noise is introduced into the system at an esti- 
mated level of 2 millimag with periodicities of 2 and ~20 
minutes. No corrections have been applied for this effect. 



2.3 OHP 1.9 m and SOPHIE 

WASP-3 was observed with the Observatoire de Haute- 
Provence's 1.93 m te lescope and the SOPHIE spectrograph 
jBouchv et alj|200fj ). over the 8 nights 2007 July 2 - 5 and 
August 27 - 30; a total of 7 usable spectra were acquired. 
SOPHIE is an environmentally stabilized spectrograph de- 
signed to give long-term stability at the level of a few ms . 



We used the instrument in its high efficiency mode, acquiring 
simultaneous star and sky spectra through separate fibres 
with a resolution of R=40000. Thorium- Argon calibration 
images were taken at the start and end of each night, and 
at 2- to 3- hourly intervals throughout the night. The radial- 
velocity drift never exceeded 2-3 m/s, even on a night-to- 
night basis. 

Conditions during both runs varied from photometric 
to cloudy, but all nights were affected by strong moonlight. 
As WASP-3 has magnitude V ~ 10.5, integrations of 900 sec 
give a peak signal-to-noise per resolution element of around 
40-50. The 2MASS colours and reduced proper motion for 
WASP-3 suggest a spectral type of about F7-8V, hence we 
cross-correlated the spectra against a G2V template pro- 
vided by the SOPHIE control and reduction software. 

In all spectra the cross-correlation functions (CCF) 
were contaminated by the strong moonlight. We corrected 
them by using the CCF from the background light's spec- 
trum (mostly the Moon) in the sky fibre. We then scaled 
both CCFs using the difference of efficiency between the two 
fibres. Finally we subtracted the corresponding CCF of the 
background light from the star fibre, and fitted the resulting 
function by a Gaussian. The parameters obtained allow us 
to compute the photon-noise uncertainty of the corrected ra- 
di al velocity measurement (any) , using the relation detailed 
in ICollier Cameron et all l|2007al '): 

(try = 1.7 * y / (FWHM)/(SNR * Contrast) 

Overall our RV measurements have an average photon-noise 
uncertainty of 14 m/s. As our radial velocity measurements 
are not photon-noise limited, we quadratically added a radial 
velocity component to those uncertainties of about 10 m/s 
(more details in Section 3.2.1). The log of the observations 
and barycentric RV is given in Table 1. 



3 RESULTS AND ANALYSES 
3.1 Stellar parameters 

The SOPHIE spectra are individually of modest signal- 
to-noise, but when summed together they are suitable for 
a preliminary photospheric analysis of WASP-3. However, 
from experience we have found that the SOPHIE standard 
pipeline reduction does not fully remove the scattered light 
component within the spectrograph. While this does not af- 
fect radial velocities significantly, it can nonetheless have 
subtle effects on absorption line depths, adversely affect- 
ing the derived spectral synthesis parameters. Therefore we 
carefully re-reduced the first three raw images taken over 
2-5 July 2007 with the reduce echelle data reduction pack- 
age (jPiskunov fc Valentill2002f ). paying careful attention to 
the issue of scattered light. These data are least affected by 
moonlight. 

F ollowing our analysis of WASP-1 (Stcmpels et alJ 
l2007f) . we employed the methodology of IValenti fc Fischer! 
(2005), using the same tools, techniques and model atmo- 
sphere grid. We us ed the iDL-based software Spectroscopy 
Made Easy (sme) l|Valenti fc Piskunovl 1 19961 ) to calculate 
and fit synthetic spectra using a multi-dimensional least 
squares approach. 

We concentrated our analysis on five regions in the spec- 
trum (see Figure [5}. These regions allow us to constrain the 



© 2002 RAS, MNRAS OOO.njfTol 



4 D. Pollacco et al. 



Table 1. Journal of radial-velocity measurements of WASP-3. The 1SWASP identifiers give the J2000 stellar coordinates of the pho- 
tometric apertures; the USNO-B1.0 number denotes the star for which the radial-velocity measurements were secured. The quoted 
uncertainties in the radial velocity errors include components due to photon noise (Section l2.3H and lOm/s of jitter (Section l3,2.1H added 
in quadrature. The fourth and fifth columns give the FWHM of the CCF dip and the contrast of the dip as a fraction of the weighted 
mean continuum level. The signal-to-noise ratio at 550nm is given in column six. 

BJD r. cxp V r FWHM Contrast S:N Notes 

(s) km s -1 km s — 1 % 



1SWASP J183431. 62+353941.4: USNO-B1.0 1256-0285133 = GSC 02636-00195 = WASP-3 
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Table 2. Parameters for WASP-3 as derived from the SME anal- 
ysis of the SOPHIE spectroscopy. 

Parameter WASP-3 



T cff 6400 ± 100 K 
logg 4.25 ±0.05 

[M/H] 0.00 ±0.20 

■usini 13.4 ± 1.5 km/s 
u ra d -5.490 km/s 



stellar effective temperature, T e g, (through the broad wings 
of Ha and, to a lesser extent, Na I D 5890A), gravity, logg, 
(through Mg I b 5175A and Na I D 5890A) and the metallic- 
ity, [M/H], (through the weak photospheric absorption lines 
m the 6000-6200A region) . We also measured the abundance 
of lithium from the Li l 6708A line. The combined spectrum 
is not of sufficient quality to perform a detailed abundance 
analysis. The parameters we obtained from this analysis are 
listed in Table [2] and a comparison between observed and 
synthesized profiles is shown in Figure [5] In addition to 
the spectrum analysis, we also used Tycho B and V, and 
2MASS photometry to estimate the effective temperatur e 
using the Infrared Flux Method l|Blackwell fc Shallis|[l977l ). 
This gives T c g = 6200 ± 300 K, which is in close agreement 
with that obtained from the spectroscopic analysis. The Ty- 
cho and 2MASS colours (V - K = 1. 32, V - H = 0.2) 
also suggest a spectral typ e of F7-8V (jStreet et al.l 120071 ; 
ICollier Cameron et al.ll2006t l. 

In our spectra the Li I A6708A line is weak, but still mea- 
surable, and we derive a Lithium abundance of logn(Li) + 
12 = 2.0 — 2.5. However, at this stellar temperature it is 
thought that stellar age does not corre late well with Lithium 
abundance |Sestito fc Randichl l2005h . so we have exam- 
ined the evolution ary tracks for l ow an d intermediate-mass 
stars presented by iGirardi et ail l|2000h using a maximum- 
likelihood fitting routine. Using our derived stellar parame- 



ters we find the stellar mass M, = 1.24 ± 0.08 with an age 
of 0.7- 3.5 Gyr. 

3.2 The reflex motion of the host star 

3.2.1 Markov-chain Monte Carlo analysis 

The SOPHIE radial- velocity data measurements are plotted 
in Figure [3] together with the best-fitting global fit to the 
SuperWASP-N, IAC80 and Keele transit photometry. Since 
the timing of the transits and the radial-velocity solution 
both provide information about the orbit, we modelled the 
transit photometry and the reflex motion of the host star 
simultaneously. 

The model of the primary star's radial-velocity orbit 
is parametrised in the usual way by the primary's radial- 
velocity amplitude Ki, the centre-of-mass velocity 7, the 
orbital eccentricity e and the longitude uj of periastron. 

The transit profile was modelled u sing the small- 
planet approximation of iMandel fc Ago] <|2002h . with the 
4-coef ficient nonlinear limb-darkening model of IClaretl 
(2000). We used i?-band limb-darkening coefficients for the 
SuperWASP-N data, whose unfiltered wavelength response 
is centred near the R band, and for the Keele .R-band 
data. We used the appropriate V and /-band limb-darkening 
coefficients for the IAC80 photometry. The transit model 
was characterised by the epoch To of mid-transit, the or- 
bital period P, the duration tr from first to fourth con- 
tact, the squared ratio AF — (R p /R,) 2 of the planet ra- 
dius R p to the stellar radius R,, and the impact parameter 
b — a(l — e cos Et) cos i/R* of the planet's trajectory across 
the face of the host star. Here a is the orbital semi-major 
axis, Et is the eccentric anomaly at the epoch of transit 
and i is the orbital inclination. The ratio of the stellar ra- 
dius to the orbital separation is then given approximately 
(or exactly for a circular orbit) by 

R* _ tr 7T 

~~a ~ T (1 + TAP) 2 - b 2 

( Collie r Cameron et al 1 l2007bf ). The orbital semi-major axis 
is derived from the orbital period and the stellar mass M* 
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Figure 2. A comparison between the observed SOPHIE spectrum of WASP-3 and the calculated spectrum obtained from spectral 
synthesis with SME, using the atmospheric parameters from Table [2] The white regions are excluded from the spectral analysis, mainly 
because of the presence of telluric absorption. Light shaded regions were used to determine the continuum level, and the remaining dark 
shaded regions to determine the stellar atmospheric parameters. All five spectral sections were used simultaneously. The five sections 
contain (a) the Mg I b triplet at 5175A (b) the Na I D doublet at 5890A/ (c) a large region with well-isolated lines of a wide range of 
metals (d) the Ho line at 6563A/ and (e) the Li 1 line at 6707A. 



via Kepler's third law. The s tellar mass is estimated from th e 
J—H colour as described by I Collier Cameron et al.l(|2007bh . 

The set of nine parameters 

{To ,P,tr, AF, b, M» ,Ki,e,u>} thus defines both the 
transit light curve and the form of the reflex velocity 
variation. We compute the photometric zero-point offset 
Am of the observed magnitudes rtlj from the model fj,j 
derived from a given set of parameters: 



Am — 



E 



3 



The weights Wj are the inverse variances l/<xf of the indi- 
vidual observations. Similarly, we compute the radial veloc- 
ity 7 of the system's centre of mass as the inverse- variance 
weighted mean offset between the observed radial velocities 
Vk and the model values for the current model parame- 



ters: 

7 = ^ • 

We quantify the goodness of fit to the data by the com- 
bined x 2 statistic for the combined photometric and radial- 
velocity data: 



£ 

3=1 



(nij - fij - Am) 2 



+£ 



7) 2 



Markov-Chain Monte-Carlo analysis has recently be- 
come established as an efficient and reliable method 



for establishing both photo metric jHolman et al ] 120061 : 
iBurke et al.ll2007h and orbital l|Fordll2006l : lGregorvll2007h pa- 
rameters of close-orbiting giant exoplanets. We determined 
the photometric and orbital parameters of the WASP-3 
system using the Markov-chain Monte-Carlo algo rithm de- 
scribed in detail bv lCollier Cameron et all l|2007bl ). to which 
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Table 3. WASP-3 system parameters and 1-<t error limits derived from MCMC analysis. 



Parameter 



Symbol 



Value 



Units 



Transit epoch (BJD) 
Orbital period 
Planet /star area ratio 
Transit duration 
Impact parameter 



T 
P 

(Rp/Rs) 2 
t T 
b 



2454143.8503 ±o' oo3 

1 8zlfi8CM+ - 000 ' 002 

l.OTUOO1_ .()()()()()2 

0106+ 0002 

U.U1UD_ Q Q004 

1H0+ 0009 
u.niu_ 001g 

-+0.051 

-0.166 



0.505^ 



days 
days 

days 
R* 



Stellar reflex velocity 
Centre-of-mass velocity 
Orbital semimajor axis 
Orbital inclination 



Ai 

7 



0.2512 
-5.4887 
0.0317 
84.4 



+0.0079 
0.0108 
+0.0013 
0.0018 
+0.0005 
-0.0010 
+2.1 
0.8 



km s 1 
km s — 1 
AU 

degrees 



Stellar mass M* 

Stellar radius i?* 

Stellar surface gravity log g t 

Stellar density p» 

Planet radius R p 

Planet mass M p 

Planetary surface gravity l°g9p 

Planet density p p 

Planet temp (A = 0) T eql 



1.24 
1.31 
4.30 
0.55 



+0.06 
0.11 
+0.05 
0.12 
+0.07 
0.03 
+0.15 
0.05 



-I ,.,+0.07 
1 - ol -0.14 

1 76+ ' 08 
i -'°-0.14 

n o 7 +0.09 

001 -0.04 

0.78t°;gf 

H33 
-76 



1960 



M 

Rq 

[cgs] 

p@ 

Rj 
Mj 
[cgs] 

P.J 
K 



we refer the reader for most details of the implementation. 
The initial photometric solution for the SuperWASP-N tran- 
sit profiles is established by our transit-search algorithm 
ijCollier Cameron et al]l200a ). The initial radial- velocity so- 
lution is an inverse variance- weighted linear least-squares fit 
assuming a circular orbit. In both cases, the initial fits also 
yield good estimates of the parameter uncertainties. The 
stellar mass is initialised at the value Mo estimated from 
the J — H colour. 

At each step in the algorithm, each of the nine proposal 
parameters is perturbed by a small random amount such 
that 

To,, = To,i_i +ff To G(0,l)/ 

and similarly for the other eight parameters. Here G(0, 1) is a 
random Gaussian deviate with mean zero and unit standard 
deviation. The adaptive step-size controller / is initially set 
to 0.5, and evolves as the calculation progresses, ensuring 
that roughly 25 percent of proposal sets are accepted. 

The prior probability distributions for all nine parame- 
ters are treated as being uniform. The parameters P, tr, 
AF, M* and K\ are required to be positive. The im- 
pact parameter and eccentricity are restricted to the ranges 
< b < 1 and < e < 1, while the longitude of perias- 
tron is restricted to the range —ir<u<n. The decision 
on whether or not to accept a set of proposal parameters is 
made via the Metropolis-Hastings algorithm using the loga- 
rithmic likelihood functional 

_ 2 (M.,i-M ) 2 (logg^-logg*) 2 

where roggr*,i is computed directly from the mass M*,, and 
radius This imposes a Gaussian prior on the stellar 



mass, and indirectly on the radius. The prior forces the stel- 
lar mass to be close to the initial estimate Mo with an as- 
sumed uncertainty au = O.IMq. The prior on loggr* en- 
sures consistency with the spectroscopically-measured value 
log g* = 4.25±0.05, and thus helps to reduce the uncertainty 
in the stellar radius if the impact parameter is not strongly 
constrained by the photometry. 

If a new set of proposal parameters yields Qt < Qi-i, 
the fit to the data is improved and the proposal is accepted. 
If Qi > Qi-i, the parameter set is accepted with probabil- 
ity exp[(Qi_i — Qi)/2\. We find that the solution converges 
within a few hundred steps to a stable, optimal solution. 
After this initial "burn-in" phase, we re-scale the photomet- 
ric error bars so that the contribution of each photometric 
dataset to % 2 ls equal to the associated number of degrees 
of freedom. For the radial velocity data, we estimate the ad- 
ditional variance needed to match their contribution to x 2 
with the number of degrees of freedom. This is equivalent to 
adding radial- velocity "jitter" with amplitude 10 ms _1 in 
quadrature with the photon noise uncertainties giving the 
values listed in Table [1] We then run the algorithm for a 
few hundred further steps and derive revised parameter un- 
certainties from the standard deviations of their respective 
Markov chains. This optimises the step length used in gener- 
ating new sets of proposal parameters. Finally the algorithm 
is allowed to run for 10 5 steps in order to map out the joint 
posterior probability distribution of the nine proposal pa- 
rameters. 

We find that for WASP-3 the nine proposal parame- 
ters sho w only weak mutual correlations. The correlation 
lengths (jTeemark et al.ll2004r i of the Markov chains for in- 
dividual parameters are typically 10 to 20 steps, so the fi- 
nal production run yielded approximately 10 4 statistically 
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Table 4. Dependence of stellar and planetary parameters on 
log g, ■ 
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independent parameter sets. In the initial runs we allowed 
all nine parameters to float, and arrived at a solution with 
eccentricity e = 0.05 ± 0.05. Since this is statistically indis- 
tinguishable from the circular orbit expected for a planet 
with such a short period, the remaining eight parameters 
were fitted assuming e = 0. 

The values of the parameters at the optimal solution 
are given, together with their associated 1-cr (68.3 percent) 
confidence intervals, in Table [3] These results are consistent 
with those derived from the spectral analysis presented ear- 
lier. 

We explored the relationship between the stellar surface 
gravity and the impact parameter of the transit by repeating 
the fit for a sequence of values of log g, with an artificially 
restricted o~i ogs = 0.01. The results are given in Table [U 
The best formal fit to the photometry is obtained for stel- 
lar surface gravities in the range 4.35 < loggr < 4.45. In this 
range, however, the stellar radius is unphysically low in rela- 
tion to the stellar mass. At the spectroscopically-determined 
loggr = 4.25, the photometric % 2 i s oni y marginally de- 
graded. At still lower values of loggr, the impact parameter 
increases to the point where the duration of transit ingress 
and egress cannot be fitted satisfactorily. Moreover, the radii 
of the star and the planet become implausibly inflated. We 
conclude that the stellar surface gravity must lie in the range 
4.25 < logg < 4.35, and the impact parameter in the range 
0.4 < b < 0.6. The limits on these parameters derived from 
the full posterior probability distribution, as listed in Ta- 
ble [3] are consistent with this conclusion. 



3.2.2 Line-bisector analysis 

It is well known that faint binaries contaminating the pho- 
tometric aperture, or even stellar activity, can influence ab- 
sorption line shape and can, in certain circumstances, mimic 
or confuse small radial velocity motions. By noting the po- 
sition of the line-bisector of the cross-correlation function, 
asymmetries in the profile will become apparent. 

We measured the asymmetries of the cross-correlation 
function peaks usi ng the line- bisector method as a function 
of rad ial velocity (|Gravl [l988). as applied bv iQueloz et al.l 
|200ll ). Figure U] demonstrates that periodic variations, in- 
dicative of line-of-sight binarity or activity are not apparent, 
and we conclude that the radial-velocity variations are gen- 
uinely due to to the orbital motion of a low mass object. 
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Figure 3. Simultaneous MCMC solution to combined 
SuperWASP-N, IAC80 V, I and Keele R photometry. The or- 
bital solution is assumed to be circular. The lower panel is the 
MCMC solution to the radial velocity data. The model here also 
shows the Rossiter-McLaughlin effect which is predicted to be 
significant in this system given the hosts vsini = 13.40 km/s. 



3.3 Adaptive Optics Imaging 

We further investigated the scenario of a triple system com- 
prising a bright single star and a faint, blended eclipsing- 
binary system by performing high-resolution H-band imag- 
ing with the near-infrared camera INGRID, fed by the 
adaptive-optics system NAOMI, on the 4.2-m William 



© 2002 RAS, MNRAS OOO.rflfTol 



8 D. Pollacco et al. 




-5.9 -5.B -5.7 -5.6 -5.5 -5.4 -5.3 -5.2 -5.1 



RV (km s"') 

08 

Figure 4. Analysis of line-bisectors in WASP-3 shows the bisec- 
tor velocity (vspan) Hoes not. correlate with stellar radial-velocity. 0.6 
This demonstrates that the cross-correlation function remains 
symmetric, and that the radial- velocity variations are not likely o.4 ; ?v 
to be caused by line-of-site binarity or stellar activity. 
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Herschel telescope. An image taken in natural seeing of 
0.8 arcsecond with corrected FWHM of 0.2 arcsecond shows 
no evidence for resolved faint companions to WASP-3. As- 
suming an F7-8V spectral type (M v ~ 3.8) for WASP-3 
would imply a distance of ~ 220 pc, hence these observations 
constrain any potential associated eclipsing binary compan- 
ion to lie within ~ 45 AU of the host. 



Figure 5. The known confirmed transiting exoplanets plot- 
ted in the Mass-Radius plane. Iso-density contours are plot- 
ted in cgs units. For clarity of scale we have not plot- 
ted HD17156 (mass 3.12 Mj, radius 1.15 Rj) or the ex- 
tremely high density object HD147506b (mass 8.04 Mj, radius 
0.98 Rj). WASP-3b is marked as a filled circle (data from 
http://obswww.unige.ch/~pont/TRANSrTS.htm and references 
therein). 



4 DISCUSSION 

In this study we have found WASP-3b to be a transit- 
ing gas-giant exoplanet with mass l.76t°'°lMj and radius 
1-311q ?4-Rj. Its host star, WASP-3, has a photospheric tem- 
perature of 6400 ± 100 K and log g = 4.25 ± 0.05, consistent 
with its F7-8V spectral type derived from 2MASS photom- 
etry. This places WASP-3b amongst the most massive of 
known transiting exoplanets (Figure [5}. Given the hosts rel- 
atively large rotational velocity and the large radius of the 
planet we would expect a significant amplitude for Rossiter- 
Mc Laughlin effect (t he m odel is included in Figure [3]) . 

ISozzetti et al.l (|2007i ) demonstrate a correlation of 
planet radius with host mass for 14 confirmed transiting 
exoplanets, and a correlation of planet mass with orbital 
period for the same sample (first noted by iMazeh et al.l 
(2005)). Figure [B] shows these relationships updated while 
in Figure [7] we also update the apparent correlation of sur- 
face g ravity with orbital period noted bv ISouthworth et al.l 
(2007). Despite the additional objects the Rp v's M, cor- 
relation remains week (even ignoring the two most massive 
objects HD17156b and HD147506b). For both the Mp v's 
P and g v's P we contend that these relationships arise 
partly through observational selection and partly through 
the effects of the intense radiation fields that these plan- 
ets are experiencing. We believe the absence of high grav- 
ity/mass bodies at longer periods is primarily a detection ef- 
fect, while the absence of low gravity /mass planets at short 
period could indeed be caused by rapid evaporation. 

The closeness of the orbit and the large radius and high 
effective temperature of the star combine to make WASP- 
3b one of the most strongly-irradiated, and hence one of 
the hottest of the known exoplanets, second only to OGLE- 
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Figure 6. Modified version of Figure 5 from ISozzetti et al.l 
(2007), incorporating a further seven newly discovered systems 
(data from http://obswww.unige.ch/~pont/TRANSITS.htm 
and references therein, but again excluding HD17156b and 
HD147506b). In each case WASP-3b is marked by the filled 
circle. 
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Figure 7. The correlation of planetary surface gravity with or- 
bital period for the 22 shor test period transiting pla netsfupdated 
version from Figure 2 from lSouthworth et all j2007l) '). 



TR-56b and comparable to OGLE-TR-132. This raises the 
possibility that the atmosphere may be hot enough for TiO 
and VO to remain in the gas phase above the temperature 
minimum, creating a hot, strongly-absorbing stratosphere 
l|Fortnev et alj|2006h . This would g ive an anomalously high 
infrared brightness temperature, as lHarrington et all (2007) 
inferred from SPITZER/IRAC secondary-eclipse photome- 
try of HD 149026b at 8/im. Being much closer and brighter 
than any of the OGLE host stars, WASP-3b is thus an ex- 
cellent candidate for future observational tests of the hot- 
stratosphere hypothesis. 
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